Stability of a Hot Smoluchowski Fluid 
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Abstract 

We study coupled nonlinear parabolic eąuations for a fluid de- 
scribed by a materiał density p and a temperaturę Q, both functions 
of space and time. In one dimension, we find some stationary Solu- 
tions corresponding to fixing the temperaturę on the boundary, with 
no-escape boundary conditions for the materiał. For the special case, 
where the temperaturę on the boundary is the same at both ends, the 
linearized eąuations for smali perturbations about a stationary solu- 
tion at uniform temperaturę and density are derived; they are subject 
to boundary conditions, Dirichlet for O and no-flow conditions for the 
materiał. 

The spectrum of the generator L of time evolution, regarded as an 
operator on L^[0, 1], is shown to be real, discrete and non-positive, even 
though L is not self-adjoint. This result is necessary for the stability 
of the stationary state, but might not be sufficient. The problem lies 
in the fact that L is not a sectorial operator, sińce its numerical rangę 
is C. 



1 Introduction 

In 16, 17, |l^, IS] we derived systems of parabolic partial differential 



eąuations as the diffusion limit of discrete stochastic models describing a 
fluid of identical particles moving in an external potential V{x). In the 
continuum limit the materiał is described by a density p{x, t) and the internal 
energy is described by a temperaturę field Q{x,t). The various models are 
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based on different kinematic assumptions, and this shows up in the details 
of the eąuations of evolution. In [|l5|, |l^ there is a maximum possible value 
of the density, called Pmaxi corresponding to the state in which all the sites 
of the discrete model are occupied. We showed that the first and second 
laws of thermodynamics hołd in these models, and that in each case, the 



system has the Onsager form |15, 2C], though it is not hnear and not near 



eąuihbrium. The difference between ||l9[| and |jT5| is that in the latter, the 
hopping probabihty in the underlying discrete model grows proportionally 
with the energy of the state the particie occupies. This assumption leads 
to a surprising phenomenon in the continuum limit: the model shows the 
Soret and Dufour effects. This had been derived in 1912 by Enskog [|l^] and 
in 1912 by Chapman [|l^] from the Boltzmann eąuation, but only for gas 
mixtures. 

These and similar models have been extended and studied by Biler, Dol- 
beault, Esteban, Hebisch, Karch, Krzywicki, Markowich and Nadzieja in a 
series of papers ||2|, ^, |, ||, |6|, |, |l^, |lj]. In these papers, V is either 
an external potential or is the mutual gravitational or electrical potential 
between the particles, and so obeys a Helmholtz eąuation with p{x, t) as its 
source. This leads to an integro-differential system, for which the authors 
have been able to prove existence and uniąueness of solutions in many cases, 
and to study the asymptotic return to eąuihbrium. 

In this paper we concentrate on the model Ę^] with the Soret and Dufour 
effects. The heat capacity per particie is unity, and so the energy density is 
naturally defined to be 

£{x,t):=p{x,t)[e{x,t)+V{x)]. (1) 

The dynamics conserves the total number of particles and the total energy, 
and this leads to the local conservation laws: 



^ + divj, = (2) 
d£ 

— + divje = 0. (3) 

The materiał current jc and the energy current je of the model are shown 
to be 

jc = -A(ev/5 + /5(i-p/p^ax)v(e + y)), (4) 

je = 2(eje-A/5(l-p/pmax)eVe) + yje. (5) 
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For rare gases, the density is much less than its maximum, and the 
system of eąuations simplifies. The currents for materiał and energy are 
then 



jc = -x{evp + pv{e + v)) (6) 
je = 2{ejc-XpevQ) + vjc. (7) 



The derivation given in starts with a Markov chain, and takes the con- 
tinuum hmit. The eąuations are therefore local in time, and if V is time- 
dependent, the same calculation shows that the right-hand side of eq. (^) 
should be pdV/dt. The model falls into the generał scheme advocated by 
Wojnarpll. 

In we show that the entropy is an increasing function of time (except 
at a stationary configuration) ; thus the model obeys both the first and the 
second law of thermodynamics if the boundary conditions allow no materiał 
or heat to enter or leave the system. The entropy is shown to be 

s = i„gp(.)d.-/(p.„„-pW)i„g(i-,/,_)d.+ 

+ j p{x)\ogQ{x)dx. (8) 
For smali /o(x)/pmax, so that we can ignore (/3/pmax)^, this simplifies to 

p{x) log p{x) dx + / p{x) \ogQ{x) dx (9) 



apart from the term / p{x) dx, which is constant in time if no materiał enters 
or leaves the system. 

Usuałly, the energy-bałance eąuation (|3|) is ignored, and the temperaturę 
put eąual to ©o for all space and time (the iso-thermal situation), and 
we reduce to the Smołuchowski eąuation for the Brownian particie in a 
potential. We avoid this simplification, which loses the Soret and Dufour 
effects that come from the cross-terms V0 in eq. (^) and 20 jc in je- The 
first is the Joule-Thompson effect (in the local form called the Soret effect), 
which ensures that a gradient in the temperaturę causes a flow of materiał, 
even when the density is constant. The term 2Qjc is the Dufour effect, which 
means that there is a flow of heat even when the temperaturę is constant, if 
there is a non-zero particie current. We might expect the flow of heat due 
to convection to be Qjc sińce the heat capacity is unity, and the factor 2 
was a surprise. In jl5U it was shown that the other hałf of the contribution 
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2Qjc is needed as the Onsager dual of the Soret term, VG, which occurs in 

jc- 

In this paper, we study these eąuations in one dimension in a bounded 
interval, which without loss of generahty may be taken to be [0, 1]. We shah 
continue to write V for the gradient, so for the rest of the paper, V = d/dx. 

Consider the the boundary conditions = 0o for all time at x = O 
and X = 1, and jc = O for all time at a; = O and x = 1. These have the 
direct physical meaning that the system is in a heat-bath at temperaturę 
Oo, and no materiał escapes through the boundary, which is impermeable to 
the particles of the fluid, but conducts heat perfectly. Such mixed Dirichlet- 
Neumann conditions for systems do not appear to have been considered 
before in that it leads to a new type of generator L of the semigroup of 
time-evolution, for which form techniąues do not work. This is because the 
operator L is not sectorial. 

It is obvious that there is a stationary solution satisfying the following 
conditions: @{x,t) = Go, and jc{x,t) = O for all space and time. It is then 
easy to solve for p in the case of a constant gravitational field (§2). We also 
briefly discuss the case when is not constant. 

In §3 we study the stability of smali perturbations about the stationary 
state, by the usual method. We limit the consideration to the case V = O, 
and so leave for the futurę the (Benard) instabilities expected for a liąuid in 
a gravitational field heated from below and cooled from above. 

We give a proof that the eigenvalues are real and non-positive, by exhibit- 
ing a 2 X 2 similarity transform (not unitary) that transforms the operator 
-L to a diagonal positive matrix, times the known negative self-adjoint oper- 
ator, the Laplacian with self-adjoint boundary conditions. This shows that 
L is closed, and leads to the generał study of operators of this type, which 
has been pursued by L. Boulton. The new difficulty is that the numerical 
rangę of the generator is the whole complex piane. 

2 Some No-flow Solutions 

Let jc and je be given by eq. (|) and (0). Let us impose the no-flow condi- 
tions dWjc = O, so in one dimension, jc is a constant. This must be zero, 
sińce we impose the boundary condition of no materiał flow at the boundary. 
Let us impose Dirichlet conditions 0(0) = 9q at x = O and 0(1) = Oi at 
X = 1. The condition for stationarity of the energy-density then leads to 
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div(p0V0) = 0; these two remarks can be expressed as 

{Qp)' + pV' = O (10) 

(ee'p)' = 0. (11) 

The second can be solved exactly for p in terms of 0. We assume that Q{x) 
and B(x)' have no zeros in [0, 1]. This gives 

P = k/{QQ'), (12) 

where k is the constant of integration, assumed not to be zero. Eliminating 
p from the first eąuation, we get 

ee" - GV' = O (13) 

This can be solved exactly for if V(x) is hnear, and = ^i- Thus 
suppose that 

V{x) = gx (14) 

as in a constant gravitational field. Then we get the obvious solution Q{x) = 
0o=const., and this leads to the solution for the density: 

p(x) = Cexp{-5x/eo}. (15) 

This is just that given by the Boltzmann factor, exhibiting the usual ex- 
ponential fall-off of the density with height. This result shows that the 
eąuations are on the right track. 

We can find many exact solutions by choosing B, satisfying the boundary 
conditions, instead of choosing V . We then just need to solve the first- 
order eąuation (p^) for V{x). Not all choices of Q{x) lead to nonsingular 
eąuations. For example, if we try to satisfy boundary conditions in which 
^0 = ^1) then a non-negative differentiable function Q[x) must have a zero 
for 0' somewhere in the interval [0, 1]. Then both p and the force — W must 
be infinite somewhere in the interval. It follows that the only continuous 
stationary solution with 6q = 9i, and 0' finite, is constant in x. It is easy 
find stationary solutions satisfying other boundary conditions; for example 
choose 0(x) = (x + 1)^; then we have a stationary solution 

p{x) = k{x + iy^ (16) 
0(x) = {x + lf (17) 
V(x) = x + x^/2. (18) 



5 



This satisfies the eąuations together with the conditions that jcix) = O for 
all X, and the boundary conditions G(0) = 1,0(1) = 4. The energy-flow, 
however is not zero, but is — 4A. This flows in at 2; = 1 and out at 2; = O, 
and is constant throughout the region. We have a driven stationary system 
far from eąuihbrium. We can follow the production of entropy exphcitly; 
heat enters the system at a high temperaturę, G = 4 and leaves at a Iow 
temperaturę, = 1. Heat leaves at the same rate as it enters. So the 
system receives entropy at the rate Q/0(1) = 4A/4 = A, and ejects it at the 
rate (5/0(0) = 4A, a net production rate of 3A. Where is the location of 
this entropy production? If we use the Onsager formahsm, we can write the 
entropy production as |l^], eąuation (4.11) 

Qs = J jc.xye + J je.xyQ. (19) 

Here, and are the thermodynamic forces, given by |jl5| (4.13) and 
(4.14). We only need 

X70 = V(l/0) (20) 
sińce jc = 0. The entropy produced by the system is therefore 

Ś= je.Xyedx = -iX[{x + l)-^]l = 3X. (21) 
Jo 

At a stationary distribution, there is no change in the microscopic probabil- 
ity distribution with time, so it seems that the von Neumann entropy (from 
which the entropy in is derived) must be independent of time. This is 
true, and is derived in eq. (4.10) in [|^]. In getting the Onsager form eq. ([l9|) , 
we discarded the boundary term, which if kept would appear as the loss of 
entropy by the system at the rate 3A, due to receiving hot heat and emitting 
cooler heat. The total change in the entropy is therefore zero, as befits a 
time-independent probabihty distribution. The environment, however, has 
an increase in entropy at the rate of 3A. 

It is easy to construct many more examples. 

3 Analysis of stability in the field-free case 

We hmit this section to the study of the case, in one dimension, on the 
interval [0,1], where V = O, and 0(0) = 0(1). We choose to write the 
dynamics in terms of the fields £ := pQ and 0; these have simple boundary 
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conditions, Neumann and Dirichlet at the points x = O and x = 1, respec- 
tively. Note that at the boundary, \/£ = O, because the materiał current, 
which is jc = —W£, vanishes on the boundary. The energy current is 
2Qjc - 2ApeVG = -2AV(G^). We can recover p = £/Q from these fields. 
The eąuations of continuity, (pi) and (H) then become 



£_ 

e 



These can be written 



e 

£ 



£& 

w 

£ 



2xe 



XV^£ 
2XV'^{e£). 



{&£)" 



£ 

2X{Q£)". 



q2 



These elear ly have the stationary solution Q{x) = G(0) = Go, £{x) 
£o, for all x € [0, 1]. To test the stabihty of this solution, we put 



e{x,t) = eo + e{x,t); 



£{x, t) = £q + e{x, t) 



(22) 
(23) 

(24) 
(25) 

m = 

(26) 



where 9{x, t) and e{x, i) are smali functions, so that we ignore powers higher 
than the first, of these and their derivatives. Up to the first order, the 
eąuations become 

e = 2A[(Go + e)(^o + e)]" 
= 2A [£oe" + Goe"] . 



Also 



2A(Go + 6) 
fo + e 



[(Go + e)(^o + e)]" 



7- ^ 

to + e 
2AQo 



[^oe" + (eo/2)e"]. 



The linearized version of the eąuations can therefore be summarised as 



2AGo 



1 7 

1/(27) 1 



L 



(27) 
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where 7 = £q/Qo. The boundary conditions are e'(0) = e'(l) = O = ^(0) = 
6{1). A necessary condition for the stability of the solution {£q,Qq) is that 
the non-zero part of the spectrum of the operator 



1 

1/(27) 



7 
1 



(8) 



(28) 



entering eq. (^) Ues in the open left-half piane. We shah prove this, and 
show that the spectrum is discrete and negative, of finite multiphcity. Zero 
is in the spectrum, sińce e{x) = non-zero constant, 0{x) = O is the corre- 
sponding eigenfunction. This perturbation corresponds to moving from one 
stationary state to another, with a different total number of particles. 

M is not symmetric, but is diagonahsable by a similarity, by Sylvester's 
criterion. In fact, there is a one-parameter family of diagonahsing matri- 
ces; for, if SMS~^ = M \s diagonal, then the one-parameter family of 
diagonal matrices B{(3) = diag[/?, commute with M, and so leave it 
diagonal: BMB~^ = M. In generał, such a transformation leaves us with 
non-self-adjoint boundary conditions. These can be expressed as a singu- 
lar perturbation of the Laplacian, in the manner of This monograph, 
however, only deals with the self-adjoint case. This extra one-parameter 
family of diagonahsing similarities can be used to ensure that the boundary 
conditions define a symmetric Dirichlet form, by a making judicious choice 
of (5. Note that it is enough to diagonalise (M — /); this leads to the special 
case, a = 1/2, of the problem of finding all matrices S such that 



and 



S 



S 



O 7 
a/7 O 



a b 
c d 



S =diag[Ai,A2 



(29) 



Direct calculation immediately gives the conditions for the off-diagonal terms 
in M to vanish: 



6^0 
ad^ 



2 2 
a 7 

2 2 
c 7 . 



To avoid the vanishing of det S we have to choose opposite signs in the 
sąuare-roots of these eąuations, for example 

b = a'ya~^^'^ 
d = —cja'^^"^. 
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Then M = diag[l + a^/^, 1 — a^^^]. We notice that this is independent of 
the value of 7. If we can show that the operator M ® dP /dx'^ is a Dirichlet 
form, it is negative sińce a = 1/2 < 1. For this to be the case, the boundary 
conditions must be self-adjoint, so that the boundary terms that arisc in the 
integration by parts vanish for functions in the domain of the operator. Let 



S 



where 



Thus, 



S = 



c —cya' 



-1/2 
-1/2 



e = ae + aja ^^'^6 



Then 



{e,9)M 



e 

9" 



The boundary terms which must vanish are thus 



+ 



+(1 - [(ce - c7a-^/2)(ce' - CTa-^/^^O]^ 
= {(1 + aV2)^a-V2a2 _ _ «i/2)^«-i/2c2| ^^^/^i _ 

For this to vanish for ah e, 9 in the domain, we need to be able to choose 
(1 + Q;^/^)a^ = (1 ~ a^/^)c^ without making S singular; this is possible. We 
conclude that our operator L is similar to a negative self-adjoint operator 
with compact resolvent, and so has non-positive discrete spectrum. The 
domain of the operator is 



{(e,^) G W^'^ : e'(0) = e'(l) = O and 0(0) = 0(1) = 0}. 



(30) 
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In order to solve the eigenvalue problem, we try 



e{x,t) = Aexp{—ujt} cos{kix) + B ex.p{—u!t} sm{kix) 
6{x,t) = C ex.p{—u)t} cos{k2x) + D exp{—Ljt} sin{k2x). 

We have proved above that co, an eigenvalue of —L, is nonnegative; ki and 
k2 must be chosen to satisfy 



UJ 



2AG(0) 



il + a'/^)kl = {l-a'/^)kl (31) 



To fit the boundary conditions, we can write (e, 9) in terms of (e, 9) and 
then use tłie given conditions. However, sińce a similarity does not alter tlie 
eigenvalue (only tlie eigenfunctions) we can avoid tłie complication of choos- 
ing S so tłiat the boundary terms vanish, and instead use any convenient 
matrix that diagonaUses M. We choose 

5=(_\ izy (32) 

This gives us the boundary conditions 

e' -9' = O, at X = O, X = 1 (33) 

e + 9 = O, at X = O, X = 1. (34) 

By imposing this condition, we get a hnear homogeneous system of eąuations 
for A, B, C, D, and the vanishing of the determinant gives us a holomorphic 
trancendental eąuation for oj. Note that the eigenvalues do not depend on 
7, the ratio of £o to ©o- We have computed the eigenvalues lj of — L, and 
indeed, they are ah non-negative. 

In it is proved that the semigroup U (t) = exp{Lt} defined by a sec- 
torial operator L with rangę in the left half-plane is a contraction semigroup 
in L^. Herę, however, the conditions for the theorem fail, and the best we 
can do in is to get the ineąuahty for u orthogonal to the one-dimensional 
subspace of stationary solutions: 

||exp{Lf}n||2 < exp{-Ait}||S'|| ||n||2 (35) 

where Ai is the smallest non-zero eigenvalue of —L. Since \\S\\ ||5'^"'^|[ > 1 in 
our model, the system is not obviously stable for smali t. This problem does 
not arise for an isolated system, in which there is no flow of heat or materiał 
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at the boundary; this gives us Neumann conditions on the boundary for both 
V£ and VG, which is a self-adjoint condition. Then the numerical rangę of 
M ® (P /dx'^ is determined by its symmetric part, which is non-positive. 
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